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Abstract. Dirichlet integrals and the associated Dirichlet statistical densities are 
widely used in various areas. Generalizations of Dirichlet integrals and Dirichlet mod- 
els to matrix-variate cases, when the matrices are real symmetric positive definite or 
hermitian positive definite, are available |4|. Real scalar variables case of the Dirichlet 
models are generalized in various directions. One such generalization of the type-2 or 
inverted Dirichlet is looked into in this article. Matrix-variate analogue, when the matri- 
ces are hermitian positive definite, are worked out along with some properties which are 
mathematically and statistically interesting. 
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1. Introduction 

This paper deals with probability densities on the space of matrices. All the matrices 
appearing in this article are p x p hermitian positive definite unless specified otherwise. 
X,Y, ... will denote matrices whose elements are functionally independent real scalar 
mathematical variables or random variables. X,Y, . . . will denote matrices whose elements 
are in the complex domain. Constant matrices will be denoted by A,B, . . . whether the 
elements are real or complex, tr(-) will denote the trace of the matrix (•), |(-)| will denote 
the determinant as well as absolute value of (•) and | det(-) | will denote the absolute value 
of the determinant of (•). Transpose will be denoted by a prime, complex conjugate by a 
bar, conjugate transpose by a star. Thus X = X* > will mean the hermitian matrix X is 
positive definite. dX will indicate the wedge product of differentials in X. For example, 
whenX = (jty) and real, 

AX = AudXij when all x, ; 's are distinct 

= f\j>jdxjj = Aj<jdxij when X = X'. 

Jo < x = x* < i(-)dX indicates the integral of (■) over all X such that X = X* > as well as 
/ — X > 0. In other words, all eigenvalues of X are in the open interval (0,1), where / 
denotes an identity manix. For any positive definite hermitian matrix A we denote by A3 
the hermitian positive definite square root of A. 
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In our discussions, we need a few Jacobians of matrix transformations and integrals 
over real scalar functions of matrix arguments. These will be stated here without proof. 
For proofs and other details, see |4|. 

Y = AXA*,\A\ ^O^d? = |det(A)| 2p dl = |det(AA*)| p dX, (1.1) 

where A is a constant matrix (p. 183 of l4l ) 

? =£ -1 , |X| ^0^dF= \det(X)\~ 2p dX forX=X*, (1.2) 

(p. 190 of E|). Let f = (tjj), Uj = for i < j, ijj = tjj > (real and positive) for j = 
1 , . . . ,p. That is, T is a lower triangular matrix with real positive diagonal elements. Then 



p 



X = ff*^dX = 2 p U] t jj jdf. (1.3) 

The complex matrix-variate gamma, denoted by T p (a), and the gamma integral are 
defined as follows: 

r p (a) = 7t Ek r n r(a)T(a - 1) • -r(a -p + 1), 91(a) >p- 1 (1.4) 
det(l)r- p e- tr( * ) dl, (1.5) 

X=X*>0 

where 9t(-) denotes the real part of (•), (p. 188 of [4|). The integral in (1.5) when evaluated 
with the help of (1.3) yields (1.4). The complex matrix-variate beta, denoted by B p (a, j3), 
and the beta integrals are defined as follows (p. 198 of |4|): 

S p (a,P) = y^iw *(") >P~ 1 ' *(P)>P-1 d.6) 

L p(U + p ) 

|det(l)| a - p |det(/-l)|' 3 - p dl (1.7) 

0<X=X*</ 

|det(l)| a - p |det(/ + l)r (a+ ' 3) dX. (1.8) 

X=X*>0 

The functions of X appearing in (1.7) and (1.8), divided by the normalizing constant 
B p (a,l5), produce the complex matrix-variate type-1 and type-2 beta densities respec- 
tively (p. 357 of 1 4 1). Dirichlet integrals and the associated Dirichlet densities come in 
naturally in order statistics problems, in reliability analysis and in certain survival anal- 
ysis problems. In Bayesian statistical analysis a Dirichlet model is usually taken as the 
prior distribution for multinomial probabilities. In random division and certain geometri- 
cal probability problems Dirichlet model comes in naturally, see for example 1 5 1 . Dirichlet 
model is extended to matrix-variate Liouville type by |2|. See [1 1 for some applications 
in engineering, [3 1 for statistical applications and |6| for some physics problems. 

The model that we are going to deal with in the present article is the complex matrix- 
variate analogue of an extended inverted or type-2 Dirichlet model incorporating succes- 
sive sums of the variables into it. 
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2. Matrix-variate analogue of an extended Dirichlet model 

Let X\,...,X k be hermitian positive definite p x p matrix random variables having the 
joint density function 

f{Xy ,...,X k ) = c k \ det(X! )\ ai ~ p ---\ dct(X k )\ 

x I det(7 + X 2 + ■ ■ ■ + X k ) |ft • • • | det(7 + X k ) |&-' 

x \fei{I + Xi + ---+X k )\- {ai+ --- +a M + ^ + ---+^ (2.1) 

for 9t(ce/) > p - 1, 3i(a j+ i H h a k+1 + J3, H h /3 fc ) > p — 1 , ; = 1,. . . ,k, and 

f{X\ ,.. .,X k ) —0 elsewhere, where c k is the normalizing constant. We will study some 
properties of (2.1) in the present article. Since 



[ ••• [f(X 1 ,...,X k )dX l A...A6X k = l, 

J X\ J X'- 



JXi JX k 

we can evaluate c k by successive integration starting with X\ . For fixed X 2 , . . . ,X k let 

U=\ Id^x^-p 

JXi=X*>0 

x |det(/ + li + ---+^)r( ai+ - +a *+ 1+ ' 3 ' + - + ^dl 1 . 

Note that 

\det(I + X l + ---+X k ) \ = \det(I + X 2 + ---+X k )\ 

x |det[/+(/ + ! 2 + "-+^r 1/2 
xX x {I + X 2 + ---+X k )- l l 2 ]\. 

Now, make the transformation 

?! = (I+X 2 + ■ ■ ■ +X k y l l 2 X x {I+X 2 + ■ ■ ■ +X k y l l 2 . 

Then 

d?! = \det(I + X 2 + ---+X k )- p \dX 1 . 
Substituting for X\ in terms of Y\ we have 

Li = |det(/ + ! 2 + "-+^)r (a2+ "' +a * +1+/3l+ "' +/y 

x / I det(y! ) I "1 -p I det(/ + ^ ) |-(«i+-+^+i+/ 3 i+-+^)dF 1 

= |det(/+X 2 + • • • +l it )|-(«2+-+^+i+/ 3 i+-+ft) 

r p (ai )f p (a 2 + ■ ■ ■ + a k +i + ft + ■ ■ ■ + ft) 
X f p (ai + --- + a fc+ i + /3i + --- + j3 fc ) 

forSR(ai) > 5R(a 2 H h a^+i +jSH \-j3 k ) > p-1. The ?i -integral is evaluated 

by using the type-2 beta integral of (1.8). Successive integrations of X 2 , . . . ,X k will yield 
the result as follows: 



-1 _ A r p(»j+i + - + %!+ A/ + ■ ■ ■ +ft) (22) 
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for%{aj)>p-\, %(a J+ i + --- + a k+1 + pj + --- + p k )>p-l, j = l,...,k. 

We can look into some interesting results and the corresponding matrix transforma- 
tions. These will be stated as theorems. We need two more results in order to establish our 
main results. These will be listed as lemmas. 

Lemma 2.1. Let X and A be p x p hermitian positive definite matrices. Then Y\ — 
A i/2£ A i/2 an dY 2 =X l / 2 AX 1 / 2 have the same eigenvalues. 

This can be easily seen by looking at the determinantal equations for the eigenvalues 
| A l/2£ A l/2 _ A/ | = q ^ |£ A _ A/ | = q ^ 1/2 _ XI\ = 0. 

Thus, Y\ and Y 2 have the same equations giving rise to the same eigenvalues X\ , . . . , X p , 
h > 0, ./' = l,...,p. 

Lemma 2.2. Let the common real eigenvalues ofY\ andY 2 of Lemma 2.1 be distinct. Then 
the wedge product dY\ — dY 2 in the integrals. 

Proof. Let U and V be unitary matrices with real diagonal elements such that 

Wi = U*YxU =D = diag(Ai, . . . ,X P ) = V*Y 2 V = W 2 . 
Then from Theorem 4.4 of |4 1 




dY { = dWi = < I 1 \X k -Xj\ z } dD AdG 




d? 2 = dW 2 = <Y\ \h ~ U } dD A dG 2 , 



and 



where dGi and dG 2 are the following: 

dGi = A[U(dU)} and dG 2 = A[V(dV)]. 

(dU) and (dV) denote the matrices of differentials (entry-wise derivatives) in U and V 
respectively. Now from Corollary 4.3.1 of |4| 

dGi = I dG 2 



ld { (p) Jo dp) F p (p) 

where <5i (p) is the orthogonal group of unitary matrices with real diagonal elements. 

Theorem 2.1. Let X\ , . . . ,X k be matrix-variate random variables having the joint distri- 
bution as in (2. 1 ). Consider the transformation 

?! =(/+*! + ■■ • +X k )- l ' 2 X l {I + X l + • • • +X k )- 1 ' 2 

Y 2 = (/ + X 2 + ■ ■ ■ +X k )- l l 2 X 2 {I + X 2 + ■ ■ ■ +X k )- l l 2 

Yk = (I + X k )- l l 2 X k {I + X k )- 1 ' 2 . (2.3) 

Then Y\,. . . ,Y k are independent, and further, Yj has a type-l beta density with the param- 
eters (ccj, 05/+H \-a k+1 + Pj-\ hj3 fe ), j=l,...,k. 
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I-Y l =I-(I + X 1 +---+X k )-^ 2 X l (I + X l + ---+X k )- 1 / 2 

= (i+x l + ---+x k )-^ 2 [(i+x l + ---+x k )-x l ] 

x (i+x l + ---+x k y^ 2 
= (/+*! + •■ ■ +x k )- l ' 2 (i + x 2 + ■■■ +x k ) 

x (I + X x + ---+X k )- l l 2 
I- Y 2 = {I+X 2 + ■ ■ ■ +X k y^ 2 (I + X 3 + ■■■ +X k ) 

x {I + X 2 + ---+X k )- l I 2 
I-Y k -i = (I+Xt-i +X k )- l ' 2 {I + X k )(I + X k _ l +X k )- 1 ' 2 

I-Y k =(I+X k )-\ (2.4) 

From the above representations of / — Y\ ,...,/ — Y k , from ( 1 . 1 ), ( 1 .2) and from Lemma 2.2, 
we can evaulate the Jacobian of the transformation in (2.3) from (2.4). 

dfi = d[(I+Xi + ■■■ +X k )- l l 2 {I + X 2 + ■■■ +X k ) 

x(/ + l 1 + ---+l it )- 1 /2] 

= d[(i+x 2 + ■■■ +x k ) l / 2 (i+x l + ■■■ +x k y l 
x(i+x 2 + ---+x k yi 2 ] 

= \det(I+X 2 + ---+X k )\ p \det{I+X l + ---+X k )\- 2p dX 1 
for fixed X 2 ,...,X k . Similarly, 

d? 2 - | det(7 +X 3 + ■■■ +X k )\P\ det(I+X 2 + ■■■ +X k )\- 2p dX 2 

and finally, 

dY k = \det(I+X k )\- 2 PdX k . 

Since the transformation in (2.3) is of a triangular nature, the Jacobian matrix will be a 
triangular block matrix with the Jacobian being the product of the determinants of the 
diagonal blocks and the Jacobian is given by, ignoring the sign, 

df l A . . . A d% = | det(7 + X { + ■ ■ ■ +X k ) |~ 2/, { det(7 +X 2 

+ ---+X k )\- p ---\det(I + X k )\-"}dX l A...AdX k . (2.5) 



246 



Joy Jacob, Sebastian George andAMMathai 



From (2.3), (2.4) and (2.5) we can compute the following product: 

k 



| J] | det(Yj)\ a i-P\ det(7 - f / )|^+i+-+«%fi+ft+-+A-p|dfi A ... AdY k 

= | 1 det(X,) 1^ J I det(7 + X 2 + • • • +X k ) | A 
x |det(/ + X 3 + • • • +l J t)| fc . . . |det(/ + l A )| A -' 

x | det(7 +Xi + • • • + |-(«i+-+°ft+i+ft+-+A)dZi A . . . A dX fc . (2.6) 

Multiplying (2.6) on both sides by we have the result since the right-hand side with 
c k is the density in (2.1) and the left-hand side with c k is the product of complex matrix- 
variate type-1 beta densities. 

It is easy to see that the converse also holds. 

Theorem 2.2 . Let the hermitian positive definite matrices Yi,...,Y k be independently dis- 
tributed as complex matrix-variate type-l beta random variables where Yj has the param- 
eters (0/,a,'+iH h (Xk+i + Pj H ^Pk)f or j= l,---,k. Consider the transformation 

in (2.3) on the space of k-tuples of hermitian positive definite matrices X\,... ,X k . Then 
Xi,...,X k have the joint density as given in (2. 1 ). 

Thus Theorems 2.1 and 2.2 also provide a characterization of the density in (2.1). It is 
known that when Yj has a complex matrix-variate type-1 beta density, then / — Yj again 
has a complex matrix-variate type-1 beta density. Thus, from Theorems 2.1 and 2.2 we 
can get two more results as corollaries. One of them will be listed here as a theorem and 
it can also be proved independently by proceeding parallel to the proof in Theorem 2.1. 

Theorem 2.3 . Let X\ , . . . ,X k have the joint density in (2. 1 ). Consider the transformation 

Zj = (7+Xj + ■ • • +X k )- l ' 2 {I+X 2 + ■■■ +X k )(I+X 1 +■■■ +X k )- l l 2 

z 2 - (7+x 2 + ■ • • +x k y l / 2 (i+x 3 + ■■■ +x k )(i+x 2 + ■■■ +x k )- 1 ' 2 

Z k = {I + X k )-\ (2.7) 

Then Zi,...,Z k are independent complex matrix-variate type-l beta random variables 
with Zj having the parameters (a ; +i H \-d k+ i+fij^ \- f5 k , (Xj),for j= l,...,k. 

Theorem 2.4 . LetX\,...,X k have the joint density in (2.1). Consider the transformation 

u l =x; l/ \i+x 2 + ---+x k )x; 1/2 
u 2 = x~ l/2 (i+x 3 + ■■■ +x k )x~ 112 



U k =X k ~\ (2.8) 

Then U\, . . . ,U k are independent complex matrix-variate type-2 beta random variables 
with Uj having the parameters (oij+i H h (X k+ \ + j3 j H h fi k , af),for j = l,...,k. 
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Proof. From eqs (2.8), (1.1), (1.2) and Lemma 2.2 we have the following: 

dUi = |det(/+X 2 + ---+X ft )Hdet(Xi)r 2p dX 1 
for fixed X2, . . . ,X k . 

dU 2 = | det(7 + X 3 + ■ ■ ■ + X k ) \p | det(l 2 ) \~ 2p dX 2 , 

and finally 

dU k =\det(X k )\- 2 PdX k . 

Since the transformation in (2.8) is of a triangular nature, we have the Jacobian given by 

dUi A . . . A dU k = | det(Xi ) \~ 2 p . . . \ dst(X k ) \- 2p \ det(/+X 2 + • • • +X k )\ p 
x\det(I+X 3 + ---+X k )\P...\det(I+X k )\P 
xdX 1 A...AdX k . (2.9) 



From (2.8), 



1/2 



/ + Oi = i + x~ 1 12 (i + x 2 + ■ ■ ■ + x k )x~ 1/2 
=x; l/2 [x 1 + (i+x 2 + ---+x k )]x; 
= x; 1/2 (i+x { + ---+x k )x; 1/2 

i+u 2 = x~ 1/2 {i+x 2 + ■■■ +x k )x~ 1/2 



I + U k =X k l/2 (I + X k )X k l/2 . (2.10) 
Now from (2.8), (2.9) and (2.10) we have 

k 



( J] \det(Uj)\ a ^ + - +a ^ l+l3 J + - +p K- p 



x | det(/ + C7,.) Jdf/j a ... A df/fc 

= |n|det(X y )|">-^| |det(/+X 2 + • • • +X k )\^ ■ ■ ■ |det(/+Xjt)|&-> 

x | det(/+Xi + • • • +l^)|-( a i+-+«*+i+ft+-+ft)dl 1 A ... A dXfc. (2.1 1) 

Multiply both sides of (2.1 1) by c k to see the result. 

The converse also holds. Thus Theorem 2.4 and its converse also provide a charac- 
terization of the density in (2.1). It is known that when Uj has a complex matrix-variate 
type-2 beta distribution then UJ 1 has a complex matrix-variate type-2 beta distribution 
with the parameters interchanged. This property also gives a couple of results. Instead of 
UJ 1 , a slightly different transformation will be considered in the next theorem. 
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Theorem 2.5. Let X\,...,X k have the joint distribution as given in (2.1). Consider the 
transformation 

v x = (i+x 2 + ■ ■ ■ +x k y l i 2 x x (i+x 2 + ■ ■ ■ +x k )- x ' 2 
v 2 = {i+x 3 + --- +x k )-^ 2 x 2 (i+x 3 + ■■■ +x k )- 1 ' 2 

%-i = {i+x k )- l l 2 x k ^{i+x k )- l i 2 
v k =x k . 

Then V\ , . . . , V k are independent complex matrix-variate type-2 beta random variables 
with Vj having the parameters (a,-, Clj+i H Vu k +\+fij-\ \-Pk)>f or j = l,---,k. 

The proof can be given by using the steps parallel to the ones in the proof of Theo- 
rem 2.1. The converse of Theorem 2.5 is also true. Further, Theorem 2.5 and its converse 
also provide a characterization for the density in (2.1). By exploiting the relationships 
between complex matrix-variate type-1 and type-2 beta random variables one can derive 
a number of results and a number of interesting matrix transformations. These will not be 
enumerated here in order to save space. 
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